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We studied the statistical properties of the luminosity distance and internal extinction data of type la supernovae in the lists 
published by Tonry et al. 1 2003 1 and Barris et al. <2004> . After selecting the luminosity distance in an empty Universe as 
a reference level we divided the sample into low z < 0.25 and high z > 0.25 parts. We further divided these subsamples 
by the median of the internal extinction. Performing sign tests using the standardized residuals between the estimated 
logarithmic luminosity distances and those of an empty universe, on the four subsamples separately, we recognized that 
the residuals were distributed symmetrically in the low redshift region, independently from the internal extinction. On 
the contrary, the low extinction part of the data of z > 0.25 clearly showed an excess of the points with respect to an 
empty Universe which was not the case in the high extinction region. This diversity pointed to an interrelation between 
the estimated luminosity distance and internal extinction. To characterize quantitatively this interrelation we introduced 
a hidden variable making use of the technics of factor analysis. After subtracting that part of the residual which was 
explained by the hidden variable we obtained luminosity distances which were already free from interrelation with internal 
extinction. Fitting the corrected luminosity distances with cosmological models we concluded that the SN la data alone 
did not exclude the possibility of the A = solution. 
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1 Introduction 

The la type supernovae are unique tools for studying several 
important large scale properties of the Universe. In particu- 
lar, their role in proving the existence of the positive value 
of the cosmological constant appeared to be fundamental. 
There are several coordinated efforts to get a large number 
of SN la of high z to give a statistically firm support for the 
nonzero cosmological constant (e.g. Riess et al. .l998..2004i 
Perlmutter et al. H^^ Tonry et al.|5003 , Barris et aL lM??! 
Astier et al. 2006; but see also Gott et al. l200U Meszaros 
120021 Rowan-Robinson f2002t . A comprehensive review on 
the cosmological implications from observations of type la 



supernovae was published by Leibundgut ( EOOli l. 

In reducing the raw observational data a crucial point is 
the estimation of the effect of the interstellar and intergalac- 
tic dust. All kinds of extinction have to be removed from the 
data before using it for testing cosmological models. The 
foreground extinction of our Milky Way is well studied and 



can be removed with certainty (Schlegel et al. I1998I I. Sev- 
eral studies indicate that the effect of the intergalactic dust 
does not play a significant role as well (Leibundgut 2001 
and the references therein). The presence of the intergalac- 



tic dust would increase the mean color excess of the distant 
objects. Unless this dust distributes very homogeneously its 
effect would increase also the scatter of the observed colors 
and magnitudes of the distant supernovae in a clear contra- 
diction to what is observed (Perlmutter et al. ll999> . The sec- 
ond property is also valid for the color-independent (grey) 
intergalactic extinction. Estimating the intrinsic extinction 
of the host galaxy poses the most serious problem. 

The observed extinction originating in the host galaxy 
depends not only on the amount and properties of the dust 
content but also on the location of supernovae within the 
host. The morphological type of host galaxies may vary on 
a wide scale from the early type ellipticals to the late type 
spirals and irregulars indicating a wide range of dust con- 
tent. There are controversial results on the average proper- 
ties of the dust in the SN la host galaxies of high redshift. 
Estimating the B-V, V-R, and V-I color excess for 20 SN la's 
Riess et al. (1996^ concluded that the ratios of selective to 
total extinction from dust in distant galaxies hosting SN la's 
are consistent with the galactic extinction law. In a recent 



paper Knop et al. J2003> did not find evidence for anoma- 
lous reddening from an independent set of 1 1 high-redshift 
supernovae observed with the Hubble Space Telescope. 
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In contrast, Clements et al ("2004) presented deep sub- 
millimetre observations of sixteen galaxies at z=0.5, se- 
lected through being hosts of a type la supernova and sug- 
gested that dust in supernova host galaxies at z=0.5 could 
produce a dimming that is comparable to the dimming at- 
tributed to accelerated expansion of the Universe. Based on 
deep submillimeter observations of 17 galaxies at z=0.5 that 
are hosts of a type la supernova Farrah et al. (2004 1 empha- 
sized the need to carefully monitor dust extinction when us- 
ing type la supernovae to measure the cosmological param- 
eters. Reindl et al. (,2005 J found that the law of interstellar 
extinction in the path length to the SN in the host galaxy is 
different from the local Galactic law. 

All these controversial results indicate that our knowl- 
edge on the properties of dust in SN la host galaxies is far 
from being complete. This situation motivated us to study 
thoroughly the statistical properties of the existing SN la 
data in particular an eventual systematical error in the avail- 
able luminosity distance and host galaxy extinction estima- 
tions. 

This paper is organized as follows: Section 2 summa- 
rizes the statistical properties of the data and claims to find 
an interrelation between the estimated luminosity distance 
and internal extinction. Section 3 deals with testing cos- 
mological models on the corrected data. Section 4 summa- 
rizes the main results of the paper The main points of this 
work were summarized previously in the paper of Hetesi 
and Balazs ( l2005t . 

2 Statistical properties of the data 

2.1 Descriptive statistics of the sample 

The most comprehensive list of SN la objects is the com- 
pilation published by Tonry et al. (|2003jl. It consists of 230 
SN la events but out of them only 188 have given internal 
host extinction. We added to this sample 23 more SN la of 
high z, published by Barris et al. (2004). The mean error 
of the estimated logarithmic luminosity distances amounted 
to 0.04 with a standard deviation of 0.02. After rejecting 
10 outliers, having uncertainties exceeding more than 3cr of 
the mean standard deviation of the estimated luminosity dis- 
tances, we got a sample size of 201 objects for the further 
statistical analysis'. The standard deviation of the errors of 
the distances did not depend on the value of extinction. Us- 
ing this fact we gave to all extinction data equal weight in 
the further analysis although their error was not given in the 
list of Tonry et al. (2003 1. The estimated host extinctions 
spread over a region of Ay = — 4.1 magnitude but except 
of a few outliers the bulk majority of data are concentrated 
in the 0-1 magnitude range. 

' Tlie compilation of Riess et al. 120041 consists of 186 SNIa objects, 
having a significant overiap with those of Tonty et al. ("2003"i. Based on this 
sample Jain and Ralston 1 2005 i pointed out a correlation between the resid- 
uals of the luminosity distances and the estimated extinction of the hosts. 
They did not attempt, however, to remove the effect of this interrelation 
from the distances. 



The distribution of the internal extinction can be eas- 
ily modelled by random orientation of a dusty galactic disc 
to the line of sight^. Assuming that the SN la supernovae 
belong to the old disc population with a scale height of 
300 pc (della Valle & Panagia 1992 Wainscoat et al. il"992l 
gave 325 pc for the old disc population) and the dust layer 
in the disc have a scale height of about 130 pc and 0.18 
mag visual extinction perpendicular to the symmetry plane 
one gets the simulated distribution as displayed in the upper 
panel of Figure[2 The lower panel of Figure[2shows the dis- 
tribution of observed extinctions for comparison. The pro- 
nounced peak at the low extinction part of the histogram is 
probable accounted for supernovae in front side of the host 
galaxy while the second less pronounced one is caused by 
those at the opposite side. The observed distribution is prob- 
able biased at the high redshift part of the sample towards 
lower extinctions because close to the limit of detection the 
heavily obscured supernovae can be missed by the obser- 
vation. Despite of the simplicity of the model applied the 
simulated data reproduce remarkably well the basic char- 
acteristics of the observed distribution within the limits of 
statistical inference. A Kolmogorov-Smirnov test revealed 
that the difference between the simulated and measured dis- 
tributions can be accounted for chance with a p — 0.547 
probability. 

The redshift distribution of the sample is bimodal (see 
Fig. There is a dip at z = 0.25 splitting the data set into 
a low and a high redshift part. In the low redshift domain 
the luminosity distances derived from various cosmological 
models differ only within the statistical uncertainty of esti- 
mating them directly from SN la events. At the same time, 
the accuracy of the measurements allows us to distinguish 
between cosmological models in the z > 0.25 range. As- 
suming that z is given, the SN la luminosity distance, de- 
rived from an empty Universe, represents an upper bound 
among those obtained from A = cosmological models. In 
what follows, we consider this bound as a reference distance 
at a given z. If the distances obtained from SN la observa- 
tion exceeded significantly this bound the introduction of a 
nonzero cosmological constant seems to be inevitable. 

Tonry et al. ( 120031) multiplied the luminosity distances 
by the present value of the Hq Hubble constant and listed 
their logarithm, along with the estimated uncertainty. 

For further statistical studies we calculated the s = {Id— 
ldo)/ai standardized deviation from an empty model, where 
Id is the measured logarithmic luminosity distance as given 
in Tonry et al. (|2p03l), ai the corresponding measurement 
error and Ido the calculated logarithmic distance in a given 
model universe. 



^ This is of course a very cmde approximation since the morphological 
types of the SN la host galaxies show a wide variety of morphologies, 
including undisturbed ellipticals, spirals and disturbed systems (Farrah et 
al. UM\ . 
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Fig. 1 Histogram of simulated SNe extinctions (top), and 
the real sample (bottom). The extinction is in magnitude 
on the horizontal axis. There is a well pronounced peak 
at Ay = due to the objects at the front side of the 
host galaxy. The second smaller peak can be accounted for 
objects at the opposite side. According to a Kolmogorov- 
Smirnov test the difference between the simulated and mea- 
sured distributions can be accounted for chance with ap = 
0.547 probability. 



-1.00 
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Fig. 2 Histogram of z distribution in the data set. The ver- 
tical dashed line indicates a cut between the low and high 
redshift part of the sample. 



For given values of J^a, i^j\/atter distances are calcu- 
lated by the following equality (Carroll et al. 1992j : 



[(1 + z'fil + nMz') -z'{2 + z')nAr'/^dz'){i) 



where f^A + + = 1 and S{x) is defined as 
sinh(x), X, and sin(a;) for ftk > 0, ilfe = 0, and il^ < 0, 
respectively. 

The s deviations of the measured distances from the 
calculated ones, based on some cosmological model, dis- 
tributed symmetrically in the case of a good fit. In the fol- 
lowing we performed several sign tests for the whole sam- 
ple and several subsamples splitting the original one with 
z — 0.25 and the median of the Ay extinction. We use the 
median instead of the sample mean value because the for- 
mer is less sensitive to outliers. Table^shows these quanti- 
ties for the two redshift range defined and the whole sample, 
respectively. 



Table 1 Statistical mean and median of the Ay internal 
extinction for the low and high z part and for the whole 
sample. Note the difference in the values between the low 
and high redshift part of the sample. The difference may be 
accounted for the observing bias in the high z part. 



redshift 


Mean 


Median 


N 


z < 0.25 


0.306 


0.250 


140 


z > 0.25 


0.194 


0.140 


61 


Total 


0.272 


0.190 


201 



2.2 Performing tests on the data 

A combination of the median of Ay internal extinction and 
the z — 0.25 cut, divided the data set into four subsam- 
ples. We used sign tests to study whether the distribution of 
the standardized residuals, introduced in the l2.1l subsection. 
was symmetric in relation to the level defined by the empty 
model, within these subsamples. A significant excess of the 
"H-" signs favors a A 7^ model and in the opposite case 
one gets support for traditional Friedmann solutions. 

We made the assumption that the probability of the "H-" 
(or "-") signs of the s residuals followed a Binomial distri- 
bution in the form of: 



(2) 



where n means the number of trials, k that of successes and 
p is the probability of one of these events (p = 0.5 in sym- 
metric case). The mean value of successes is given by np. 
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One may compute the probability of the \np—k\ < 6 devia- 
tion from the mean assuming that it happens only by chance; 



P{\np-k\ <d)^ 



E 

-k\<5 = 



Pi 



(3) 



The formula given by Equation (|3 has the advantage 
that it is exact and works correctly even in the case of a 
small sample size. We carried out sign tests on the subsam- 
ples defined by cutting the data set at the median in the low 
and high redshift part, separately. Tables |2] and |3l summa- 
rize the results. As Table|2]demonstrates the subsamples are 
symmetric, independently of the internal extinction. On the 
contrary, the results summarized in Table |3 indicates a dif- 
ference in the statistical properties of the low and high ex- 
tinction part of the subsample (see Figure|3- While the low 
extinction part clearly shows an excess of the points above 
the reference line it is not the case among the points above 
the Ay median. This fact indicates some sort of interrelation 
between the internal extinction and the standardized residu- 
als in the high z domain. 



Table 2 Sign tests for the low redshift (z < 0.25) part of 
the sample, divided by the median value of the Ay internal 
extinction. The last column (Pi) gives the probability that 
the numbers of "-" and "H-" signs differ only by chance. The 
distribution of the standardized residuals are symmetric, in- 
dependently on the value of Ay. 





k ("-") 


n-k C+n 


n 


Pi 


Av < 0.25 


34 


35 


69 


1.000 


Av > 0.25 


34 


37 


71 


0.813 


Total 


68 


72 


140 


0.800 



Table 3 Sign tests for the high redshift (z > 0.25) part of 
the sample, divided by the median value of the Ay internal 
extinction. The last column (Pi) gives the probability that 
the numbers of "-" and "H-" signs differ only by chance. Un- 
like to Table |2]the low extinction part of the sample has an 
excess of the "H-" signs which can not be accounted only by 
chance. On the contrary, the distribution in the Ay > 0.140 
domain is still symmetric. The difference in the behavior be- 
tween the Ay < 0.140 and Ay > 0.140 part of the sample 
indicates an interrelation between the s normalized residu- 
als and the Ay internal extinction. 





k ("-") 


n-k ("+") 


n 


Pi 


Av < 0.140 


9 


21 


30 


0.043 


Av > 0.140 


17 


14 


31 


0.720 


Total 


26 


35 


61 


0.305 



The sign test is appropriate for studying the significance 
of the excess of points above the s — line and could give 
a clear support for A 7^ cosmology. One may wonder, 
however, whether additional statistical tests also confirm the 



Fig. 3 The extinction - standardized residual relation of 
the low (top) and high-z (bottom) subset. Vertical dashed 
lines mark the median of the extinction data. Horizontal 
dashed lines mark the reference level of an empty Universe. 
Note the difference between the left and the right panel. In 
the z < 0.25 case the distribution of residuals is symmetric 
to the reference level of an empty Universe, independently 
from the extinction. On the contrary, the low extinction part 
(left from the median line) of the z > 0.25 panel clearly has 
an excess of the points above the reference line but it dis- 
appears at higher extinction values displaying a pronounced 
negative trend. 



result obtained from the sign test. To proceed in this way 
we performed three additional tests on the high z part of 
the sample: Student's t, Mann-Whittney and Kolmogorov- 
Smirnov tests. The Student's test compares the means while 
the others the distribution of s at both side of the median of 
Ay. We summarized the results in Table|2 
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Table 4 Additional tests on the high z part of the sample where <tav ; cTs are the standard deviations of Ay and s, re- 
spectively. For each data point the value of the / variable 
(the factor score) is resulted with keeping the first PC and 
EA ; £s are the residuals when substituted / into Equation 
0. 



Type of test 


sign. 


No. of objects 


Student's t 


0.044 


61 


Mann-Whittney 


0.025 


61 


Kolmogorov-Smirnov 


0.005 


61 



One may infer from Table |4] that the results of the ad- 
ditional tests confirm the conclusion obtained from the sign 
test. 

2.3 Pearson's correlation and factor analysis 

To quantify the interrelation pointed out in subsection 12.21 
we computed Pearson's linear correlation between Ay and 
s in the low and high z domain, separately. Inspecting Table 
|5]infers that there is a very significant correlation in the high 
redshift part unlike the low redshift one. 

We made the assumption there was an / hidden variable 
which was present in both s and Ay and it was responsible 
for the interrelation. From strictly statistical point of view, 
it was not necessary to specify the physical nature of this 
variable. Following this assumption, the observed s and Ay 
could be expressed in terms of / by the following system of 
equations: 



Ay 
s 



Ao 
So 



/ + 



(4) 



where Aq, sq are constants and ea, £s represent noise terms. 
To estimate Aq, sq and / we invoked the standard tech- 
nique of factor analysis. In many statistical software pack- 
ages (e.g. in the widely used SPSS^) the default solution of 
a factor model is equivalent to the computation of principal 
components (PC) of the correlation matrix of the observed 
variables {Ay and s in our case). For obtaining PCs one has 
to solve the eigenvalue problem of the correlation matrix of 
the observed variables: 



1 r 
r 1 



= A 



(5) 



where r is the Pearson correlation between Ay and s, 
{01,12} are the components of the eigenvector, and A is 
the eigenvalue of the matrix. 

If one used principal components for establishing a fac- 
tor model the usual way is to make a distinction between the 
'significant' and 'non-significant' PCs and keeping only the 
former ones. A common procedure for making this kind of 
distinction is to keep only those eigenvectors which belong 
to eigenvalues above a given threshold. The widely used 
Kaiser criterion (Kaiser [T960t keeps eigenvalues of A > 1. 
We used this criterion in our analysis. 

The {fli , a2 } components can be used to obtain the con- 
stants in Equation according to the following equalities: 



^0 — (^AvO-l 



So = 0-502 



(6) 



SPSS is a registered trade mark for Statistical Package for Social Sci- 



Table 5 Pearson's linear correlation between the internal 
extinction and the s standardized residual. The third column 
gives the probability for the purely random occurrence of 
the correlation. 





CORRiAv,s) 


Sig. (2-tailed) 


N 


z < 0.25 


-0.068 


0.427 


140 


z > 0.25 


-0.411 


0.001 


61 


Total 


-0.122 


0.086 


201 



The / background variable reproduced well both the Ay 
internal extinction and the s residual as displayed in Figure 
|4] One has to remove the effect of this variable from the 
luminosity distances before using them for testing cosmo- 
logical models. 



2.3.1 Caveats 

The high level of significance obtained for the Pearson's 
correlation has to be treated with some caution. Usually, 
one assumes that both variables have Gaussian distribution 
when calculating the significance level of the Pearson's cor- 
relation. In our case this assumption is far from the reality in 
particular at Ay in the z > 0.25 range. There are 57 objects 
in the < Ay < 0.5 part and only 4 with 0.5 < Ay < 1.0. 

If one considered only the data points in the < Ay < 
0.5 range of the z > 0.25 sample the Pearson's correlation 
is still negative but drops back to -0.221 with a significance 
of 0.094. 

To get a measure of correlation independent on the as- 
sumption of normality we split the sample into four parts 
with the medians of s and Ay. The number of data points in 
these quadrants form a 2 x 2 contingency table. Making use 
this table one can compute other measures of correlation as 
listed in Table |6l 

Table 6 Different types of correlation in the high redshift 
part of the sample 



Type of corr. 


CORR 


Sig. 


N 


Pearson's linear (0 < Ay < 1.0) 


-0.411 


0.001 


61 


Pearson's linear (0 < Ay < 0.5) 


-0.221 


0.094 


57 


Pearson's R (2 x 2 cont. table) 


-0.381 


0.005 


57 


Spearman (2x2 cont. table) 


-0.381 


0.005 


57 



One may infer from Table|6lthat measure of correlation 
obtained from the contingency table is close to the Pearson's 
linear correlation computed for the whole z > 0.25 sample. 
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Fig. 4 Dependence of the Ay internal extinction (upper 
panel) and the s residual (lower panel) on the / background 
variable obtained from the factor analysis. Full lines mark 
the relationships between the observed and hidden variables 
as obtained from the factor analysis and dashed ones de- 
fine the 95% confidence region for the fit. The effect of the 
hidden variable has to be removed from the luminosity dis- 
tances before using them for testing cosmological models. 

3 Testing cosmological models 

As we pointed out at the end of Section|2]the luminosity dis- 
tances had to be freed from the interrelation with the internal 
extinction, before testing cosmological models. The factor 
analysis splits the actual residual into two parts: s = s' + s" 
where s' represents that part explained by the / background 
variable while s" stands for the true one which is indepen- 
dent from the amount of the internal extinction. 

Going back to the definition of s as given in Subsection 
12.11 we may write Id — Ido + ai{s' + s"). Subtracting s' 
from this equation, that part explained by the background 
variable, we obtain: Id'^ = Id — cris' = Mq + ais" which is 
already free from the interrelation with internal extinction. 
This luminosity distance can be used for testing cosmologi- 
cal models. 



The luminosity distance of an object of z redshift can be 
calculated in the framework of a given cosmological model. 
Besides z it depends on f^A and fl^j cosmological param- 
eters. It is widely accepted, based on the recent results of 
SN la data and the WMAP project, that f^A = 0.7 and 
— 0.3 favoring a f^A + = 1 flat Euclidean model 
(Hansen et al. l2004l Verde l200^ Wright 2003 1. 

We did not intend to repeat the estimation of JIa and 
riM using the original data of Tonry et al. J2003> and Har- 



ris et al. J2004> . To make a comparison we used only our 
corrected data. For this purpose we calculated the dis- 
tance between our corrected data and those obtained from 
cosmological models according to the equation below: 



X 



iid--id-'''%z,nA,QM)y 



1=1 



(7) 



Minimizing this distance with respect to JIa, gives the 
best estimate of these parameters according to the observed 
set of Id'^. Assuming that the W — Id'^"-^'^ residuals follow a 
Gaussian distribution and introducing the L = log P likeli- 
hood function we obtain: 



L = - 



constant 



(8) 



The confidence of the parameters can be estimated making 
use of the following theorem : 



2(L, 



^o) « xl 



(9) 



where L^ax , Lq mean the likelihood function at the max- 
imum and the true value of the parameters estimated, re- 
spectively (for the details see Kendall & Stuart ll973t . The 
degree of freedom is p = 2 in our case. 

Based on our corrected luminosity distances we calcu- 
lated the best fitting values of Ha, by minimizing in 
Equation and, consequently, maximizing L in Equation 
(|8ji at the same time. 

As a best fit we obtained JIa = 0.47 and Hm — 0.43 
which considerably differ from the canonical f^A = 0.7 and 
riM — 0.3. In Figure|5lwe displayed the scatterplot between 
z and the deviation of the observed luminosity distances and 
the calculated ones in an empty Universe before (left panel) 
and after correction (right panel). It is worth noticing the 
much smaller scatter of the corrected data around the best 
fitting line in the z > 0.25 region in comparison with the 
non-corrected data. Computing the value for one degree 
of freedom we obtained Xdf — 1-46 for the non-corrected 
sample and Xdf — 0-38 for the corrected one. 

Equation (|9jl enabled us to calculate the confidence in- 
terval of the parameters estimated. After fixing a value for 
Xp Equation ^ specifies a boundary in the {^Ia, ^m} pa- 
rameter space. The particular value specified for Xp gives a 
/3 probability that the true value of the estimated parameters 
are outside and I — (3 that inside this boundary. 1 ~ (3 gives 
the confidence level of this region. Following this procedure 
we calculated the confidence and the result is shown in the 
right panel of Figure|6l In the left panel of this figure we also 
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Fig. 5 Deviation of the observed logarithmic SN la lumi- 
nosity distances from those in an empty Universe as a func- 
tion of the z redshift. Upper panel shows this relationship in 
the case of uncorrected data and the lower one after correc- 
tion. The best fitting solutions are also indicated. Note the 
much smaller scatter of the corrected data around the best 
fitting line in the z > 0.25 region in comparison with the 
non-corrected data. Blue and red lines mark the best fit to 
the corrected (ilj^ = 0.47 ; J7m = 0.43) and the canonical 
model (i^A = 0.7 ; ^Im = 0.3), respectively. Green line 
indicates the A = Euclidean model. 

give the confidence of the parameters using the uncorrected 
sample for comparison. 

A comparison of the left and the right panel of Fig- 
ure |6l clearly shows the basic difference between the pa- 
rameters obtained from the uncorrected and corrected data. 
While the uncorrected data strongly supports the A / 
solution A = is well within the 95% confidence region 
in the corrected case. It is also important to note that the 
canonical iljy ~ 0.7 and 57m = 0.3 solution is outside 
of the 99% confidence region of the corrected data. How- 
ever, the 51a + — 1 line referring to an Euclidean Uni- 
verse crosses the la confidence region with a best estimate 
of 51a = 0.55 and 51 a/ — 0.45 values. 

The statistical results obtained in Section |3 are purely 
phenomenological and give no hint of the probable reason 
of the interrelation between the luminosity distance and in- 
ternal extinction as given in the data set of Tonry et al. 
(12003 1 and Barris et al. (2004 1. A detailed careful simula- 
tion of the procedure of obtaining the luminosity distance 
and internal extinction would be desirable to uncover the 




Fig. 6 Confidence region of 51a, in case of the orig- 
inal data (left panel) and the same after correction (right 
panel). The outer boundaries of the black, orange and blue 
areas represent 67%, 95% and 99.3% confidences, respec- 
tively. Red line separates the accelerating and decelerating 
and green line the open and closed model. Note that unlike 
the uncorrected case, the corrected data are consistent with 
the existence of a A = solution. 

probable reason of the interrelation found in Section |2] A 
detailed study of this kind may reveal whether this interrela- 
tion has some astrophysical reason or is purely a byproduct 
of the procedure for obtaining the luminosity distance and 
internal extinction. 

4 Summary and Conclusions 

We studied the basic statistical properties of the SN la sam- 
ple published by Tonry et al. (I2003t and Barris et al. ( I2004t . 
The observed distribution of the internal extinction due to 
the host galaxies can be well modelled by a dusty disc ori- 
ented randomly to the line of sight. We divided the sample 
into the low and high z part at a dip in the distribution at 
z = 0.25. 

We further divided the low and the high redshift part of 
the sample by the median of the internal extinction. Four 
subsamples resulted in this way. We selected the redshift 
- luminosity distance relationship of an empty Universe as 
a reference level and calculated the s standardized devia- 
tion of the logarithmic luminosity distances of Tonry et al. 
( I2003> and Barris et al. ((20041 from those of the empty 
model. 

Performing sign tests on the standardized residuals, on 
the four subsamples separately, revealed that s is distributed 
symmetrically in the low redshift (z < 0.25) part of the 
sample, independently on the internal extinction. On the 
other hand, the low extinction part of the data of z > 0.25 
clearly showed an excess of the points of s > which was 
not the case among those of having extinction above the me- 
dian. This diversity in the behavior of the points below and 
above the median pointed to an interrelation between the s 
residual and the Ay internal extinction. 

To characterize quantitatively the interrelation between 
s and Av we computed the Pearson's linear correlation in 
the low and the high redshift part, separately. Assuming that 
the interrelation can be represented by a hidden common 
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variable in s and Ay we constructed a factor model for 
explaining the observed variables. We used PCA and the 
Kaiser criterion for verifying the factor model. 

After subtracting that part of the residual which was ex- 
plained by the factor model we obtained a corrected sample, 
free akeady from the interrelation between s and Ay - Min- 
imizing the distance between the corrected data and a 
cosmological model with respect to ft^, Hm one may con- 
clude that there was no need for A ^ 0. 

We have to emphasize, however, than one came to this 
conclusion assigning the same weight to the low and higher 
extinction part of the sample in the z > 0.25 region. Nev- 
ertheless, some concerns can be made for the reliability of 
the estimated extinction in the host galaxy. Without the pro- 
posed correction the low extinction part of the z > 0.25 
sample clearly supports the A 7^ cosmology. 

Our result was purely phenomenological based only on 
the statistical properties of the data. Consequently, further 
detailed studies are required whether the interrelation be- 
tween the luminosity distance and the internal extinction in 
the high z part of the sample has an astrophysical mean- 
ing or is simply a byproduct of the way as Id and Ay were 
obtained from the data. 
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